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ABSTRACT 



We present photometry and analysis of the microlensing alert MACHO 96-LMC-2 
(event LMC-14 in Alcock et al. 2000rj| ). This event was initially detected by the 



MACHO Alert System, and subsequently monitored by the Global Microlensing Alert 
Network (GMAN). The ~ 3% photometry provided by the GMAN follow-up effort 
reveals a periodic modulation in the lightcurve. We attribute this to binarity of the 
lensed source. Microlensing fits to a rotating binary source magnified by a single lens 
converge on two minima, separated by A% 2 ~ 1. The most significant fit XI predicts a 
primary which contributes ~ 100% of the light, a dark secondary, and an orbital period 
(T) of ~ 9.2 days. The second fit X2 yields a binary source with two stars of roughly 
equal mass and luminosity, and T = 21.2 days. 

Observations made with the Hubble Space Telescope (HST)Q resolve stellar neighbors 
which contribute to the MACHO object's baseline brightness. The actual lensed object 
appears to lie on the upper LMC main sequence. We estimate the mass of the primary 
component of the binary system, M ~ 2 M Q . This helps to determine the physical size 
of the orbiting system, and allows a measurement of the lens proper motion. For the 
preferred model XI, we explore the range of dark companions by assuming O.IMq and 
1.4Mq objects in models Xla and Xlb, respectively. We find lens velocities projected to 
the LMC in these models of vxia = 18.3±3.1 kms -1 and vxu = 188±32kms _1 . In both 
these likelihood analysis suggests an LMC lens is preferred over a Galactic halo 

lens, although only marginally so in model Xlb. We also find vx2 = 39.6 ± 6.1 kms -1 , 
where the likelihood for the lens location is strongly dominated by the LMC disk. In 
all cases, the lens mass is consistent with that of an M-dwarf. Additional spectra of 
the lensed source system are necessary to further constrain and/or refine the derived 
properties of the lensing object. 

The LMC self-lensing rate contributed by 96-LMC-2 is consistent with model self- 
lensing rates. Thus, even if the lens is in the LMC disk, it does not rule out the possibility 
of Galactic halo microlenses altogether. Finally, we emphasize the unique capability of 
follow-up spectroscopic observations of known microlensed LMC stars, combined with 
the non-detection of binary source effects, to locate lenses in the Galactic halo. 

Subject headings: dark matter - gravitational lensing - stars: low-mass, brown dwarfs - 
binaries: general 



lr The NASA/ESA Hubble Space Telescope is operated by AURA, Inc., under NASA contract NAS5-26555. 
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1. Introduction 

The interpretation of gravitational microlensing results towards the Magellanic Clouds (e.g., 



Afonso et al. 1999; Alcock et al. 2000b) has been hindered by the unknown location of the lensing 



systems. Exceptions to this are caustic crossing binary lens events MACHO LMC-9 ( [Alcock et al 



2000a ), where a sparsely resolved caustic crossing suggests an LMC lens, and MACHO 98-SMC-l 



(Alcock et al. 1999a; Afonso et al. 2000| ), where a lens association with the SMC is more certain. 



In this paper we present an additional "exotic" microlensing event seen towards the Magellanic 
Clouds, MACHO 96-LMC-2. 

The lightcurve of MACHO 96-LMC-2 exhibits deviations from the standard microlensing fit 
similar to those that are expected if a binary source is lensed ( |Griest fc Hu 1992j ). This effect is the 



"inverse" of the parallax effect due to motion of the Earth around the Sun flGould 1992 ; Alcock 



et al. 1995), and may be referred to as the "xallarap" effect, where the orbital motion occurs at 
the lensed source. Detection of the xallarap modulation in a microlensing lightcurve allows us to 
fit the semi-major axis of the orbiting system in units of the lens' projected Einstein ring radius. 
An estimate of the physical semi-major axis of the system then allows a second constraint (along 
with the event timescale t) on the 3 degenerate lens parameters mass, distance, and transverse 
velocity. [Han fc Gould (1997|) describe the use of this effect in discriminating between Galactic 
halo and LMC lenses. 

A detection of this type of modulation in a microlensing lightcurve is possible within a certain 
range of event parameters. First, the binary source should have an orbital period similar to or 
shorter than the event timescale, such that the sources accelerate appreciably during the time 
they are microlensed and an orbital period may be determined. Second, the orbital separation of 
the sources should not be much smaller than the lens' Einstein ring radius projected to the source, 
otherwise the system appears to the lens as essentially a single object. This biases the detection of 
the xallarap effect towards events where the lens is close to the sources, analogous to the parallax 
effect, which is most easily detected when the lens is relatively close to the Sun-Earth system. A 
comprehensive review of microlensing with rotating binaries (lenses, sources, and observers) may 
be found in Dominik (1998; ), 



2. Observations 

Microlensing Alert MACHO 96-LMC-2 was detected and announced on Oct 3, 1996, with the 
MACHO object at an observed magnification of A ~ 1.8. The source for this event is located 
at a = 05:34:44.437, 5 = -70:25:07.37 (J2000), in the south-east extreme of the LMC bar. This 
object was constant at V = 19.42 ± 0.15, (V - R) = -0.03 ± 0.10 in ~ 700 observations over the 
4.2 years preceding this brightening, where these magnitude errors are dominated by uncertainty 
in the calibration of the MACHO database ([Alcock et al. 1999b| ). The MACHO ID number for 



this star is 11.8871.2108. A 25" x 25" postage stamp, taken from the MACHO R-band template 
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observation of this field and centered around the lensed object, is presented in Figure |T|. 

Nightly observations were requested on the CTIO0 0.9m telescope as part of the Global 
Microlensing Alert Network (GMAN) microlensing follow-up program QAlcock et al. 1997b| ; |Alcock| 



et al. 1999a} [Alcock et al. 2000a| ). Data were obtained in both B and Kron-Cousins R for the 



duration of this event. Final sets of baseline observations were made ~ 800 days (~ 8 times the 
duration of the event) after the peak. An additional set of observations was made at the UTSO 
0.6m telescope. However, closure of the telescope in 1997 prevented baseline measurements from 
this site. 

After the event, cycle-7 HST images were taken with the target star centered in the Planetary 
Camera of HST instrument WFPC2. These integrations included four 500 second exposures in 
each of three bands V, R, and I. The images were combined using the IRAF routine imcombine, 
along with a sigma clipping algorithm to remove cosmic rays. Aperture photometry was performed 
on all stars using the DaoPhot package (Stetson 1994a), with centroids derived from a PSF fit. 



We used a 0.25" aperture, and corrected to a 0.5" aperture using the brightest stars in the field. 



We corrected for the charge transfer effect and calibrate the magnitudes using the Holtzman et al 



(1995) calibrations. A portion of the combined R-band image of the WFPC2-imaged field around 



event 96-LMC-2 is presented in Figure |lj. 

The MACHO/GMAN data for 96-LMC-2 are presented in Figure |. The MACHO data were 
reduced with MACHO 's standard photometry package SoDOPHOT, with minimum errors of 
0.014 added in quadrature. The CTIO and UTSO data were simultaneously reduced with the 
ALLFRAME package ([Stetson 1994bj ), and the error estimates are multiplied by a factor of 1.5 to 
account for global systematics (such as flat-fielding errors and the swapping of CCD detectors) in 
the time series of data. 



3. Microlensing Fits 

3.1. Standard Microlensing 

Standard microlensing fit parameters, including the effects of unlensed contributions to the 
source baseline flux, are presented in Table [I]. These include 

• to, the time of closest approach of lens to source-observer line of sight, 

• t, the Einstein ring diameter crossing time, 

• u mim the lens impact parameter in units of the Einstein ring radius. 



2 Cerro Tololo Interamerican Observatory, National Optical Astronomy Observatories, operated by the Association 
of Universities for Research in Astronomy, Inc. under cooperative agreement with the National Science Foundation. 
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We also include 5 additional blending parameters, one / for each passband of observations. 
In each case, / represents the fraction of the object's baseline flux which was lensed. The x 2 /d.o.f 
for this fit is 0.91, formally an acceptable fit. However, there are periodic residuals around this 
smooth fit, especially in the CTIO data. We have plotted these residuals for the CTIO R and B 
passbands in Figure and Figure H, respectively. 



3.2. Binary Source (Xallarap) Microlensing 

Fitting this event to an orbiting binary source does provide significant improvement over the 



standard blended fit. We follow the formalism of Dominik (1998] ) for the binary source solution. 
In particular, we use the fit parameters: 

• tb, the time of closest approach of the lens to the source system center of mass, 

• t~E, the lens' Einstein radius crossing time, 

• b, the lens' impact parameter with respect to the source system center of mass, in units of 
the lens' Einstein radius projected into the source plane, 

• a, the angle between the lens trajectory and the x source axis, 

• /i, the total binary flux fraction of source 1, 

• rn-i, the total binary mass fraction of source 1, 

• p, the orbital semi-major axis in units of the lens' projected Einstein radius, 

• (3, the orbital inclination, 

• T, the orbital period in days, 

• £oj the orbital phase at time tb, 

and one / for each passband of observations, an additional 8 parameters compared to the 
standard blended microlensing fit. Finally, we assume zero eccentricity circular orbits for the 



sources, meaning inclination angle 7 from pominik (1998| ) is redundant, and is set to zero in these 
fits. 

We find 2 minima in this parameter space, separated by A% 2 ~ 1- Our most significant model 
is labeled XI, and our second most significant is X2. These fits are a further Ax 2 = 72 from 
the standard microlensing fit, which we are extremely unlikely to arrive at by chance, even given 
our additional 7 constraints. Xallarap fit parameters are presented in Table ||. Fit XI indicates 
a primary which contributes ~ 100% of the light, a dark secondary, and an orbital period of 
T = 9.22 ± 0.21 days. The second fit X2 yields a binary source of similar mass and brightness 
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stars, and T = 21.2 ± 0.54 days. These fits and the residuals of these fits around the standard 
microlensing fit are plotted along with the data for the CTIO R and B passbands in Figure || and 
Figure respectively. 

We have investigated whether or not the sources contribute significantly different flux 
fractions in the red and blue, by allowing different f\ fractions for the red and blue passbands. 
Providing this additional constraint leads to a A% 2 = — 0.24(— 0.50) for fit XI (X2), indicating this 
improvement is formally significant at only the 37% (52%) confidence level. For the fit XI class, the 
secondary source is dark in both red and blue, to within the reported accuracy of our photometry. 
For the X2 class, the best fit f\ r j f\b is 1.02. In the following we assume f\ r = fib = f\. 

The binary source fits for this event are plotted with the data in Figure Due to the different 
blend fractions for these fits, the observed object magnification (as opposed to the lensed source 
magnification) is plotted as a function of time. It is apparent that without the GMAN follow-up 
photometry, there would be little support for the binary source interpretation. In fact, 80% of the 
A% 2 between the standard fit and fit XI is contributed by the GMAN data. 



3.2.1. Colors of the Lensed Objects 

The source object's brightness is well constrained with our V, R, and I HST images. We use 
the image subtraction method of [Tomancy &: Crotts (199(f ) to locate the lensed source in the 
MACHO images to within 0.1" (2 PC pixels). Image registration allows us to uniquely determine 
the centroid of the lensed source in the HST image. This lensed flux aligns with an object with 
V = 19.46 ± 0.02, (V - R) = 0.00 ± 0.03, (V — I) = 0.12 ± 0.03, where the errors represent the 
quadratic sum of Poisson noise and an adopted error of 2% for our HST magnitudes. Within 1" 
of this source there are at least 2 neighbors, which contribute about 20% of the flux within this 
region. 

Color-magnitude diagrams (CMDs) incorporating ~1800 objects from the Planetary Camera 
chip of the WFPC2, and a region surrounding the lensed object in the MACHO focal plane, are 
displayed in Figure || The lensed object identified in each respective CMD is indicated with the 
filled circle. We have corrected the CMDs for reddening, using a characteristic LMC reddening in 
the bar of E(B - V) = 0.07 ± 0.01 (e.g., |Qlsen 1999| and Holtzman et al. 1999|) and the Landolt 



extinction coefficients of Schlegel, Finkbeiner, &; Davis (1998). The intrinsic magnitude and colors 



of the source star are V = 19.23 ± 0.04, (V - R) = -0.04 ± 0.03, (V — I) = 0.02 ± 0.03. 

The CMDs in Figure || indicate the source lies very close to the main sequence in this region 
of the LMC. There appears to be a slight excess of flux in the I passband, as can be seen in the 
(V — R), (R — I) and (V — I),V diagrams. There is no apparent excess in the (V — R), V diagram, 
suggesting this is an infrared excess. Given the direction of the reddening vector in the V, R, and 
/ passbands, this marginal excess cannot be a feature of reddening. This could possibly be due to 
the lens itself, or in the context of model XI, a signature of the dark companion to the primary. 
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We also indicate in Figure | the region 0.75 mag below the lensed object with an open circle, which 
would be the location of a single component of this object if it were a blend of equal brightness 
stars, as model X2 suggests. 



3.2.2. Properties of the Binary System 

We next attempt to estimate the mass of an object with the colors determined above. Since 
this object appears on the upper main sequence, it is likely to have higher metallicity than the 
majority of the field stars in this location. Here we consider Girardi et al. (2000) isochrones for 



z = 0.008. For model class XI, the best fits to the colors and apparent magnitude of the star 
come from objects with log(age[yr]) = 8.7 — 8.8, M = 2.1 ± O.IMq. In the case of model X2, we 
require a single object 0.75 mag dimmer than our observed object. Using the same isochrones, 
we find a wider range in acceptable star age, log(age[yr]) = 8.5 — 8.8, and a range in mass of 
M = 1.9 ±0.1M Q . 

The original fits to model XI yielded rh\ ~ f\ ~ 0. In the case of f\ = 0, we are not able 
to fit model parameters fh\ and p individually, but can only constrain their product rh\ * p. 
Therefore the fitting process was re-run, setting fh\ =0.1 {fh\ = or rh\ = 1 imply non-existent 
secondaries) and fitting for p. Therefore fit XI has 1 more degree of freedom than X2. The product 
rh\ * p was similar to that found in the original model. To explore the range of companions to 
the 2.1Mq primary, we consider a "light" dwarf secondary of O.IMq, henceforth fit Xla, and a 
"heavy" white dwarf or neutron star secondary of 1.4Mq, fit Xlb. Making these assumptions fixes 
fh\ and allows us to extract an associated p. The resulting parameters fh\ and p for fits Xla and 
Xlb are also listed in Table |2[ 

Knowing the total mass of the system and orbital period allows us to solve for the physical 
Keplerian parameters of the binary source. We find semi-major axes of 0.11, 0.13, and 0.23 AU, 
and circular velocities of 132, 154, and 120 kms -1 , for fits Xla, Xlb, and X2, respectively. 



3.2.3. Constraints on the Lensing Object 

Xallarap fit parameter p relates the scale of the lens's projected Einstein radius (-Re) to the 
binary's semi-major axis. The derivation above allows us to express this value in AU, and the 
results are presented in Table ||. This effectively leads to a measurement of the lens proper motion 

v± Re m 

To find the velocity of the lens projected to the LMC (v), we assume a LMC distance modulus 
of 18.5, or a distance of 50 kpc. We further assume the source is 1 scale height behind the midplane 
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of the LMC disk (see below) . The values of \i and v for each of our fits are presented in Table |3[ 
The error bars on these values do not incorporate uncertainty in the distance to the LMC. 

The proper motion measurement allows a one-parameter family of solutions relating the lens 
mass and distance 



4G 1 — x 

where x is the ratio of lens to source distances, D[/D s . 

With some knowledge of the kinematics of sources and lenses, we can estimate a most likely 
distance for the lensing object, given our measured v. This directly leads to an estimate of the 
lens mass from Equation (|2|). This is performed in a manner similar to Alcock et al. (1997b| ) and 



Alcock et al. (2000a) for each of the results from models Xla, Xlb, and X2. We use the following 



representation of the Galactic-LMC system: the Galactic halo is represented by Model S of Alcock 
et al. (1997a1 ); we use the preferred LMC disk model of pyuk, Dalai, &: Griest (2000[ ) , with a 
scale-height of 300 pc, tilt of 30 deg from face-on, central surface density of 190 Mo/pc 2 , and 
we assume the source to be 1 scale height behind the midplane; the Solar velocity in Galactic X, 
Y, Z is (9, 231, 16) kms -1 ; the LMC velocity is (53, -160, 162) kms" 1 , with isotropic random 
velocities of 22 kms -1 (e.g., Graff et al. 2000] ) in each direction. For the LMC halo, we assume a 



similar profile to the Galactic halo, with a central density of 0.0223Mqpc -3 , core radius of 2kpc, 
and velocity dispersion of 50kms _1 . We truncate this halo at llkpc. We note the details of the 
Galactic and LMC models do not qualitatively alter our conclusions. 

We find in each analysis that a lens residing in the LMC is preferred to a Galactic halo lens, 
although only marginally so in fit Xlb. In particular, for fit Xla the likelihood of measuring our 
value of vxia = 18.3 kms -1 is dominated by the LMC disk model. For this fit, a LMC disk lens 
is 7 times more likely than a LMC halo lens, and > 950 times more likely than a Galactic halo 
lens. We also point out that secondaries less massive than O.IMq become increasingly unlikely, as 
decreasing the mass of the secondary leads to a lower lens v. In this model, vxia already is drawn 
from the low velocity tail of the LMC disk probability distribution. As an example, a secondary 
of approximately Jupiter mass would imply a lens velocity projected to the LMC of ~ 0.2 kms" 1 , 
clearly in contradiction to the LMC disk likelihood profile in Figure §. For fit Xlb, we find 
vxu = 188 km s -1 , which is most likely to come from our model of the LMC halo, although it is 
only 4 times as likely as a Galactic halo lens. In this case, a LMC disk lens is ruled out at high 
confidence. We note that such a LMC halo population has yet to be detected directly. Given the 
broad range in secondary mass, and hence the binary's semi-major axis, explored between fits 
Xla and Xlb, we conclude this model favors a lens associated with the LMC. With our model 
X2, we find vx2 = 39.6 ± 6.1 kms -1 , whose likelihood is strongly dominated by the LMC disk. 
Figure |6| shows the probability distributions of v for each component of our Galactic-LMC model. 
Our measured v for each fit is also displayed with 1 — a errors. Only in model Xlb does a Galactic 
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halo lens appear reasonable. 

Our likelihood analysis also yields a probable distance to the lens, and from Equation (Q), a 
mass. These parameters are listed in Table || for each model. Model Xla allows the lightest lens 
mass, 0.057 +°;°3gM , while models Xlb and X2 imply heavier lenses, 0.24+J;Jjj and 0.28 ^j^M©, 



respectively. We have referenced the low mass isochrones of Girardi et al. (20001 ) to determine the 



expected brightness of a main sequence M-dwarf lens. We use the isochrones for z = 0.004, and 
an assumed lens age of 4 Gyr. We find for objects of M = O.3M , and also for the upper 1 — a 
confidence limit M = O.6M , absolute V magnitudes of 9.9 and 7.6, respectively. At the distances 
implied by the likelihood analysis, even the most luminous configuration provides negligible flux 
from the lens compared to the apparent brightness of the source object. A O.6M lens at 45kpc 
leads to an apparent lens magnitude of V = 25.9, or 0.2% of the source brightness. The apparent 
I-band brightness of this lens is / = 24.7, or less than 0.1% of the source brightness in this band. 
The secondary source is therefore a more likely origin for the possible I-band excess seen in 
Figure |5[ The small lens proper motion of ^ 1 milli-arcsecond year -1 implies it will take the next 
generation of space telescopes to be able to resolve and image the lensing object. 

We note that the lens' location (and therefore its mass) can be more accurately determined 
with direct spectral observations of the source. In particular, the superposition of primary and 
secondary spectra can provide a discriminant between models XI and X2, and lead to a better 
mass estimate for each object. Similarly, the orbital period can be precisely determined with 
radial velocity measurements. These observations will constrain the Keplerian parameters more 
directly than our (necessarily) more complicated microlensing analysis has allowed. 



3.3. Alternate Models 

In order to gauge the robustness of these conclusions, we have considered several alternate 
models, including placing the sources at larger distances behind the LMC, and increasing the 
amount of reddening to the sources. 

I 3 

Zhao (2000) predicts a strong excess in lensed source reddening for models where the LMC 
is strongly self-lensing. In these models, the source stars must preferentially lie at the back 
side of the LMC, and thus will have high interstellar reddenings. On the other hand, if the 
lenses are Galactic dark matter, the reddenings of the source stars will be statistically the same 
as surrounding stars along the same lines of sight. Alcock et al. (2000c) use HST colors of 8 



microlensed source stars to rule out, at the 85% confidence level, a model where all sources are 
located ~7 kpc behind the LMC disk. However, it is possible that roughly half of the sources 
may be in the background. Therefore, we consider a model where the lensed sources are 7 — 9 kpc 
behind the LMC. This leads to an increase in source mass of ~ O.1M , and results in projected 
velocities of vxia = 18 kms -1 , vxu = 184kms _1 , and vx2 = 40kms _1 . The v likelihood profiles 
for this source location differ significantly from those seen in Figure @. Most notably, the LMC 
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disk and LMC halo profiles allow a similar range in v, with the LMC halo generally preferred over 
the LMC disk, and median v ~ 100 km s -1 . For the lowest projected velocity (model Xla), the 
lens is consistent with neither the LMC nor the Galactic likelihood profiles. The LMC halo is the 
most likely location of the lens in all models, followed by the LMC disk in model X2, and by the 
Galactic halo in model Xlb. 

It is also possible the reddening to the sources is larger than our adopted value of 
E(B — V) = 0.07. For example, the foreground reddening maps of Schwering &; Israel ( 199l] ) 



indicate E(B — V) = 0.12, while Zaritsky (1999 ) finds typical LMC source star reddenings of 



E{B — V) ~ 0.10. If we assume a higher reddening to the source object of E{B — V) = 0.12, its 
intrinsic brightness and colors are V = 19.06, (V — R) = —0.08, (V — I) = —0.05. For model XI, 
the primary source would weigh ~ 2.3M Q , and for model X2 both sources would weigh ~ 2.1M Q . 
Here we find vxia = 17kms _1 , vxib = 181kms _1 , and vx2 = 41kms _1 , whose relative likelihoods 
can be evaluated using Figure ||. This change has more important implications for the source mass 
than for the projected velocities, and the likelihood results are similar to those from the standard 
analysis. 



4. Conclusions 

We have measured and characterized a periodic modulation in the lightcurve of microlensing 
event MACHO 96-LMC-2. We model this event using a single object microlensing a rotating 
binary source, which provides considerable improvement over a single source model. Possible 
alternate explanations to this modulation include a single source lensed by a rapidly rotating 
binary lens, or some variant of stellar pulsation. We do not consider such models. MACHO 
96-LMC-2 is not the only time that binary source effects have been detected in a microlensing 
event. EROS-II event GSA2 ( Derue et al. 199"s| ) exhibits a similar modulation around the standard 



microlensing fit. Derue et al. (1999|) also find a binary source model degeneracy similar to that 



between our models XI (dominant source) and X2 (equal brightness sources). 

We are able to constrain the projected Einstein radius of the lens (-Re) m the most significant 
of our fits to be between 0.54 and 5.5 AU, and its velocity projected to the LMC v to be between 
18.3 and 188 kms -1 . The weakness of this constraint (an order of magnitude!) is due to our 
lack of knowledge of the mass of the secondary component of the lensed binary system. In this 
model we have no means to constrain the secondary's orbit due to its negligible contribution to 
the system's brightness. We chose example secondaries separated by an order of magnitude in 
mass (O.IMq to 1.4M Q ), leading to order of magnitude constraints on and v. However, in both 
cases, a LMC lens is preferred to a Galactic halo lens. For the larger velocity lens (model Xlb) 
the object should come from an as yet undetected LMC halo population, and a Galactic halo lens 
is not strongly ruled out. The lack of direct evidence for an LMC halo population indicates model 
Xlb is best able to constrain the location of the lens to be out of the LMC disk. Our second most 
significant model leads to Re = 1.24±0.18 AU, and v = 39.6 ± 6.1 kms -1 . This model also prefers 
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a lens in the LMC, with the likelihood heavily weighted towards the LMC disk. All 3 of these 
models suggest a sub-solar mass lens, consistent with an M-dwarf star. Our derived values of Re, 
v, M[, and D\ are in good agreement with the characteristic properties of LMC lenses presented 
in Table 1 of [Han & Gould (1997Q . 



Alternate models for the source system, including one where the sources lie 7 — 9 kpc behind 
the LMC, and one where they are reddened by an additional AE(B — V) = 0.05, were also 
considered. By placing the sources far behind the LMC, we are unable to discriminate between 
the LMC disk and LMC halo population of lenses, and can only state that the lens is most likely 
associated with the LMC system. Additionally, the qualitative conclusions about the location of 
the lens are found not to be overly sensitive to the amount of reddening to the sources. 

The identification of MACHO 96-LMC-2 with a LMC lens population is consistent 
with the expected LMC self-lensing signal. This single event microlensing optical depth is 
t a = 1.1 x 10" 8 ,r B = 8.5 x 10~ 9 for criteria "A" and "B" in |Alcock et al. (2000b| ), respectively. 



This is to be contrasted with an expected self-lensing optical depth of r = 1.6 x 10 for the LMC 
disk and r = 1.7 x 10~ 8 for the MACHO component of our model LMC halo ( [Alcock et al. 2000b[ ). 



It is also possible that MACHO LMC-9 is due to LMC self-lensing ( [Alcock et al. 2000a| ), although 
we caution this interpretation is based on a caustic crossing resolved with only 2 observations, and 



other interpretations are possible. LMC-9 is excluded from event set A in Alcock et al. (2000b ), 
and has an optical depth of tb = 9.3 x 10~ 9 . The expected LMC self-lensing signal due to these 2 
events is thus likely to lie within the range 1.1 x 10~~ 8 < r < 1.8 x 10 -8 . This is approximately half 
of the expected LMC self-lensing rate. Thus the combined optical depths alone do not implicate 



the LMC as the host for the majority of the microlenses, as originally suggested by Sahu (1994 ). 

I I 

Kerins & Evans (1999) take the ensemble of events LMC-9 and 98-SMC-l and argue that 

the existence of even a few LMC self-lensing events suggests an LMC halo interpretation of all 
LMC microlensing. In this context, the appearance of another apparent LMC self-lensing event 
strengthens the case for significant LMC self-lensing. However, because of the importance of our 
result for the interpretation of LMC microlensing, we emphasize the bias a study such as this has 
against revealing a lens residing in our Galactic halo. In particular, a lensed LMC binary source is 
preferentially more likely to show xallarap modulations if the lens is also in the LMC. To identify 
Galactic halo lenses, a spectroscopic study of all Magellanic Cloud microlensed sources should be 
undertaken to assess whether or not the lensed object is, in fact, a binary system. In the case of a 
positive detection, we may set limits on xallarap modulation in the microlensing lightcurve, and 
from this a lower limit on the proper motion of the lensing object. As can be seen from Figure ^, 
Galactic halo lenses start to dominate the likelihood near v = 230 km s -1 . A lower limit on v 
near this value would be highly suggestive of a true halo lens, and thus detection of at least one 
component of the Galactic dark matter. We suggest such a study of the high magnification event 
MACHO 99-LMC-2 ( Pecker 1999|) , whose lightcurve is well covered by the MACHO/GMAN, 
MOA, MPS, and OGLE microlensing teams. 
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Fig. 1. — Each image represents a 25" x 25" field centered on the lensed source object. The image 
on the left is from the 300 sec MACHO R-band template observation, and on the right is a 4 x 500 
sec combined HST R-band image. The lensed source is in the center of each image, and is indicated 
by the circle. North is up and east is to the left. 
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MACHO Event 96-LMC-2 
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Fig. 2. — The light curve of event 96-LMC-2. The panels show the observed brightening as a 
function of time, with passbands and sites as indicated. The solid line represents the best "xallarap" 
fit XI, and the dotted line fit X2. Insets are provided to better view the region with the strongest 
binary source signal. The UTSO data are only plotted with the XI fit curve. 
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Fig. 3. — The CTIO R-passband light curve of 96-LMC-2. The top panel displays the 3 fits 
performed on the complete time-series of data, along with the CTIO R-band data in 1-day bins. 
The dotted line represents the best standard fit, and the solid and dashed lines the most significant 
and second most significant xallarap fits XI and X2, respectively. The bottom panel indicates the 
residuals of each fit and the data around the standard microlensing fit. 
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Fig. 4. — Same as Figure 0, except with the CTIO B-band light curve. 
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Fig. 5. — Color-magnitude and color-color diagrams of the field surrounding event 96-LMC-2. 
The lensed object is indicated with a filled circle. The HST CMDs are constructed from cycle-7 
observations with the target star centered on the Planetary Camera chip of the WFPC2. The 
MACHO CMD is taken from a subsection of the template observation of this field. Error bars 
indicate the uncertainties in each respective CMD. This object resides near the upper main sequence 
of the LMC field. We indicate with an open circle the region of the CMD where 2 equal brightness 
stars must be drawn from to create the observed object. 
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Velocity projected to LMC. 




Fig. 6. — Probability distributions of v for each component of our Galactic-LMC model, described 
in Sec. 3.2.3| . The measured v for models Xla, Xlb, and X2 is also indicated, with respective I — a 
errors. Note in each case, a LMC lens is preferred over a Galactic halo lens, although this preference 
is weakest for model Xlb. 
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Table 1. Standard Microlensing Fit Parameters 



Fit Parameter 


X 2 1 d.o.f. 


1872.31 / 2060 




1767.51 (7) 


t 


98.7 (40) 


^min 


0.301 (17) 


JmACHOr 


1.00 (13) 


ImACHOb 


1.00 (13) 


fcTIO R 


0.99 (12) 


fcTIO B 


1.02 (13) 


fuTSO R 


0.76 (33) 



a (JD - 2448623.50). 

Standard point source, point 
lens microlensing fit parameters 
for event 96-LMC-2. For each 
passband, / represents the fraction 
of the object's baseline brightness 
which was lensed. Reported 
uncertainties in the final significant 
digit (s) are the maximum extent 
of the surface in parameter space 
which has a % 2 greater than the 
best-fit value by 1. 
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Table 2. Xallarap Microlensing Fit Parameters 







.A. X cx YY _L U 


X9 


X 2 1 d.o.f. 


1799.9 / 2054 




1800.9 / 2053 


h a 


1767.600 (79) 




1767.80 (14) 


2 * 


101.8 (42) 




108.8 (53) 


b 


0.287 (16) 




0.246 (19) 


a 


-2.15 (37) 




0.236 (77) 


fi 


0.000 (29) 




0.52 (19) 


rh\ 


0.1 b 


0.045 0.40 


0.53 (19) 


P 


0.095 (15) b 


0.208 (33) 0.0237 (39) 


0.188 (26) 


P 


1.00 (30) 




1.43 (14) 


f 


9.22 (21) 




21.22 (53) 


lo 


-0.23 (37) 




0.40 (12) 


fMACHO R 


0.94 (9) 




0.80 (8) 


fMACHO B 


0.95 (9) 




0.80 (8) 


fcTIO R 


0.93 (9) 




0.79 (8) 


fcTIO B 


0.96 (9) 




0.81 (8) 


fuTSO R 


0.58 (18) 




0.86 (46) 



a (JD - 2448623.50). 

b For fit XI, we are only able to constrain the product of rh\ and p. 
We assume 0.1M Q and 1.4M Q dark companions to the 2.1Mq primary 
to determine these parameters for fits Xla and Xlb, respectively. The 
sources for fit X2 are estimated to be 1.9M . See Sec. |3.2.2 for further 
details. 

Xallarap microlensing fit parameters for event 96-LMC-2. The 
parameters are as defined in Sec. 3^. For each passband, / represents 
the fraction of the objects baseline brightness which was lensed. 
Reported uncertainties in the final significant digit (s) are the maximum 
extent of the surface in parameter space which has a x 2 greater than 
the best-fit value by 1. 
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Table 3. Properties of the Lensing Object 



Xla Xlb X2 



Re (AU) 


0.535 (88) 


5.50 (91) 


1.24 (18) 


\x ( kms -1 kpc -1 ) 


0.364 (62) 


3.73 (64) 


0.79 (12) 


v (kms -1 ) 


18.3 (31) 


188 (32) 


39.6 (61) 


Mi (Mq) 


5.7+™ x 10- 2 


2.4+^ x 10" 1 


2.8+f-J x 10" 1 


D t (kpc) 


49.61°;* 


no r +6.5 

38.5 _ 15? 


49.63J 



Characteristics of the lensing object, determined from our xallarap fit 
parameters and estimates of the individual masses of the lensed binary 
system. i?E represents the lens' Einstein ring radius projected to the 
source system, and v represents the lens velocity projected to the source 
system. Lens properties M[ and Di are estimated using a maximum 
likelihood technique described in Sec. 3.2.3 , and a LMC distance of 50 
kpc. 



